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Resonance  Scattering  of  Lyman- Alpha  Radiation  by 

Hydrogen  in  the  Ground  State 


1.  INTRODUCTION 

The  calculation  of  the  resonance-scattering  cross  section  for  Lyman-alpha 
radiation  is  important  for  a  variety  of  reasons.  First  of  all,  such  scattering  is 
extremely  important  in  astrophysics,  since  hydrogen  is  the  most  abundant  element 
in  "empty"  space.  Hydrogen  is  also  important  in  the  higher  portions  of  the  earth's 
atmosphere,  where  it  acts  as  a  scatterer  of  solar  radiation. 

Another  reason  for  calculating  this  cross  section  is  that  this  is  the  simplest 
resonance  cross  section  which  one  can  calculate,  since  the  wave  functions  are 
simple  and  since  a  two-level  approximation,  involving  the  ground  state  and  the 
degenerate  first  excited  state,  is  apt  to  be  especially  good.  One  may  expect  some 
insight  into  the  methods  of  calculating  resonance  cross  sections  for  more  compli¬ 
cated  atoms. 

A  more  fundamental  reason  is  that  the  resonance  line  shape,  which  is  identical 
to  the  line  shape  of  Lyman-alpha  radiation  emitted  by  excited  hydrogen,  lead 
directly  to  the  self-energy  problem  of  quantum  electrodynamics.  As  is  well- 
known,  the  position  of  resonance  (or  maximum  intensity  in  the  case  of  emission) 
is  not  precisely  at  the  Lyman-alpha  frequency,  but  is  shifted.  In  the  traditional 
calculations  of  this  shift,  the  dipole  approximation  to  the  electromagnetic  matrix 


(Received  for  publication  5  December  1973) 
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elements  is  used,  and  the  shift  is  infinite.  Since  the  advent  of  renormalization 
theory  one  uses  renormalization  to  obtain  finite  results.1 2 3 4 5 6  Where  resonance  scat- 
tering  is  treated  directly,  it  is  found  that  the  shift  is  of  the  order  of  the  Lamb 
shift  of  the  energy  levels.  However,  to  obtain  a  finite  result,  renormalization 
theory  is  used  which  includes  a  higher  order  correction. 

3 

In  the  present  paper,  we  use  Dirac's  theory  of  resonance  scattering.  In  this 

theory  only  the  matrix  elements  which  have  a  resonant  denominator  in  the  usual 

perturbation  theory  expansions  are  assumed  to  be  non-zero.  One  is  then  led  to  a 

Hamiltonian  whose  eigenfunctions  can  be  obtained  exactly.  The  Dirac  resonance 

theory  is  studied  from  a  rigorous  mathematical  point  of  view  by  K.O.  Friedrichs, 

4 

and  we  are  greatly  influenced  by  his  paper. 

However,  we  must  modify  the  Dirac  resonance  theory  somewhat  to  take  into 
account  the  degeneracy  of  the  first  excited  state  of  hydrogen.  In  our  treatment  we 
shall  consider  only  the  matrix  elements  of  the  electromagnetic  interaction 

Hj  =  i  (eft/ Me)  A  •  V  .  (1) 

where  M  is  the  electron  mass,  which  connects  the  following  states  with  each  other 

3 

by  the  emission  of  a  photon  from  the  vacuum:  IS,  2S,  2  P.  In  contrast  to  the 
traditional  treatments,  we  shall  not  use  the  dipole  approximation  but  instead  we 

5 

shall  use  exact  matrix  elements.  It  will  be  seen  that  in  contrast  to  earlier  treat¬ 
ments,  the  shift  in  position  of  resonance  will  be  finite  without  renormalization. 
Furthermore,  this  shift,  when  translated  into  shift  '  in  energy  level,  is  close  to 
the  Lamb  shift  for  the  ground  state.  We  include  all  matrix  elements  of  the  two- 
level  system,  those  which  lead  to  "forbidden"  (in  the  dipole  approximation)  as  well 
as  those  which  lead  to  "permitted"  transitions.  When  we  first  carried  out  the  cal- 

g 

culations  for  the  frequency  shift  using  the  techniques  of  the  present  paper,  we 
ignored  the  "forbidden"  transitions.  This  calculation  led  to  a  frequency  shift 


1.  Stenholm,  S.  (1973)  Quantum  theorj  of  electromagnetic  fields  interacting  with 

atoms  and  molecules,  Physics  Reports  6C(No.  1).  North-Holland  Publishing 
Co. ,  Amsterdam.  ** 

2.  Low,  F.  (1952)  Natural  line  shape,  Phvs,  Rev.  88:53. 

3.  Dirac,  P.  A.  M.  (1947)  The  Principles  of  Quantum  Mechanics,  3rd  Edition, 

Clarendon  Press,  Oxford,  p. 201. 

4.  Friedrichs,  K.O.  (1948)  On  the  perturbation  of  continuous  spectra,  Commun. 

Appl.  Math.  1:361. 

5.  Moses,  H.E.  ( 1 97 3)  Photon  wave  functions  and  the  exact  electromagnetic 

matrix  elements  for  hydrogenic  atoms,  Phys.  Rev.  A,  8:1710, 

6.  Moses,  H.E.  (19G6)  Resonance  scattering  of  Lyman-a  radiation  by  hydrogen 

in  the  ground  state.  Tech.  Note  1966-14,  Lincoln  Lab. ,  MIT. 
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somewhat  greater  than  one-half  of  the  Lamb  shift.  [The  line  shape  of  emitted 

Lyman-alpha  radiation  was  also  calculated  in  7  ,  also  taking  into  account  only  the 

permitted  transition  between  the  IS  state  and  the  23 9P  states,  using  the  exact 

matrix  elements  (calculated  differently  but  agreeing  with  other  matrix  elements  '  ) 

but  with  an  unconventional  treatment  of  the  electromagnetic  interaction.  Exactly 

6 

the  same  shift  as  found  previously  was  obtained.  ]  The  inclusion  of  matrix  ele- 

3 

ments  corresponding  to  "forbidden"  transitions  between  the  2S  and  2  P  states 

O 

and  between  the  2  P  states  themselves  gives  a  substantial  contribution  to  the 

3 

shift.  On  the  other  hand,  only  the  permitted  transition  from  the  IS  to  the  2  P 
state  contributes  to  the  natural  line  breadth.  It  would  be  interesting  to  see  what 
contributions  the  n=3  states  and  two-photon  intermediate  states  give.  There 
should  be  no  difficulty  in  calculating  these  contributions  other  than  "book-keeping.  " 
If  the  contributions  are  small,  it  would  appear  that  the  two-level,  single-photon 
theory  is  indeed  a  good  approximation. 

In  short,  the  calculations  of  the  present  paper  are  of  the  kind  that  could  have 
been  done  in  the  early  1930's  had  the  exact  matrix  elements  been  known.  The 
shift  which  we  obtain  corresponds  to  the  raising  of  the  ground  state  by  an  amount 
of  6297  MH  in  frequency  terms.  The  Lamb  shift  calculated  using  renormalization 

O 

theory  is  8126  MH.  Thus  our  calculation,  which  is  completely  free  of  infinities, 
gives  in  excess  of  75  percent  of  the  Lamb  shift.  Since  retardation  gives  such  a 
large  proportion  of  the  Lamb  shift  which  is  usually  attributed  to  renormalization, 
it  would  appear  that  renormalization  theory  should  be  re-examined  for  bound  states 

Q 

at  least.  It  ought  to  be  mentioned  that  the  inclusion  of  the  effects  of  retardation  is 
being  treated  in  recent  times  by  a  number  of  writers  from  a  Green's  functions 
point  of  view  (see,  for  example,  10).  However,  the  objectives  and  methods  are 
very  different  from  those  of  the  present  paper. 


2.  THE  DIFFERENTIAL  AND  TOTAL  CROSS-SECTIONS. 

THE  RADIATION  PATTERN  OF  THE  SCATTERED  RADIATION. 
COMPARISON  WITH  PHENOMENOLOGICAL  THEORY 


The  quantity  k  =  (2ir/X)  is  the  wave  number  of  the  incident  and  scattered  radia¬ 
tion  which  is  assumed  to  be  near  the  wave  number  n  of  Lyman-alpha  radiation. 


7.  Stroud,  C.  R. ,  Jr.  (1970)  Quantum  and  Semiclassical  Radiation  Theories. 

University  Microfilms,  Ann  Arbor. 

8.  Bethe,  H.A.,  and  Salpeter,  E.  E.  (1957)  Quantum  Mechanics  of  One-  and 

Two-Electron  Atoms,  Academic  Press,  Inc. ,  Hew  York,  p,  103. 

9.  Hoffman,  H. S. ,  and  Moses,  H.E.  (1972)  The  ultraviolet  convergence  of  the 

ground  state  of  hydrogen,  Letters  al  Nuovo  Cimento  4:54. 

10.  Gavrila,  M. ,  and  Costescu,  A.  (1970)  Retardation  in  the  elastic  scattering 
of  photons  by  atomic  hydrogen,  Phys,  Rev.  A  2:1752. 
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From  the  Bohr  formula  for  the  energy  levels  of  hydrogen 


a  =  |  a  (2) 

where  a  is  the  first  Bohr  radius  of  the  hydrogen  atom  and  or  is  the  fine  structure 
constant.  The  polarization  of  the  radiation  is  most  conveniently  described  by  the 
helicity  variable  P=±  1 .  If  0=  1 .  the  radiation  is  circularly  polarized  in  the  direc¬ 
tion  opposite  to  that  of  propagation,  whereas  if  /3=  -1 ,  the  radiation  is  circularly 
polarized  in  the  direction  of  propagation. 

We  shall  assume  that  incident  radiation  of  wave  number  k  and  helicity  /? 
propagates  in  the  positive  direction  along  the  z-axis  from  minus  infinity  and 
strikes  the  hydrogen  atom  in  the  ground  state  located  at  the  origin.  Let  us  denote 
by  T)  the  unit  vector  in  the  direction  of  propagation  of  the  scattered  radiation,  and 
let  us  introduce  the  usual  polar  coordinates  0  and  d>  by 

r)  =  (sind  cos  6,  sind  sin  6,  cos  d).  (3) 


We  shall  denote  the  cross  section  for  the  scattering  of  the  photon  with  the 
same  circular  polarization  as  the  incident  radiation  in  the  direction  2  by  Oj(d). 
This  and  the  other  cross-sections  are  independent  of  0.  Then,  Oj  (d)dfl  (where 
dQ  =  sin  d  dd  d0)  gives  the  relative  intensity  of  the  radiation  scattered  in  the  solid 
angled  ft.  Explicitly, 

cr.(d)  =  ~2~  F(k)  cos4  |  ,  (4) 

1  4  fT  * 

where  F(k)  is  the  resonance  function  given  by 


F(k)  = 


(k-K  -6l2+  y2 


and  in  which  the  half-breadth  y  and  the  resonance  shift  6  have  the  values 


6  =  -  1.  320 


y  =  -1.  046  X  10 


■2  -1 
cm 


In  terms  of  frequency  these  quantities  are 


||  =  -6297  MH,  ||  =  -49. 9  Mil. 


(5b) 
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The  differential  cross  section  for  scattering  of  radiation  whose  circular 
polarization  is  opposite  to  that  of  the  incident  radiation  is  denoted  by  anc* 

is  given  by 

a  Ad)  =  — F(k)  sin4 -I .  (6) 

*  4fC i  1 

The  differential  cross  section  for  scattering  of  ohotons  of  either  polarization  but 
with  the  incident  circular  polarization  still  given  by  /3  is  denoted  by  a^,(0)  and 
is  just  the  sum  of  cr  j (f3 )  amd  ^(O): 

<rT(fl)  =  0^0)  +  cr2(0). 

One  easily  verifies  that 

aT(0)  =  —^2  F(k)  [l  +  cos20J  .  (7) 

8k 

If  the  incident  radiation  is  unpolarized,  one  averages  over  both  values  of  (3  to 

obtain  the  cross  section  for  the  scattering  of  unpolarized  incident  polarization 

without  distinguishing  the  polarization  of  the  scattered  radiation.  Since  cT(0)  is 

independent  of  0,  a^(6)  also  gives  the  differential  scattering  cross  section  for 

unpolarized  \ncident  radiation.  It  is  of  interest  to  note  that  the  0 -dependence  of 

OrpW)  is  identical  to  that  predicted  by  the  classical  theory  of  resonance  fluores- 
1  11 

cence.  In  the  quantum  theory  of  scattering,  we  shall  see  that  this  dependence  is 
due  to  the  fact  that  in  our  approximation,  only  photons  with  angular  momentum 
quantum  number  j=  1  are  scattered. 

The  total  cross  section  E  is  given  by 

E  =  j*  aT<0)dft 

where  the  integration  is  taken  over  the  entire  sphere.  We  have 

E  =  F(k)  .  (8) 

We  can  confront  Eq.  (8)  with  the  phenomenological  results  of  Mitchell  and 
1 2  12 

Zemansky.  It  is  assumed  in  **  that  in  resonance  scattering  a  photon  is  absorbed 

11.  Panofsky,  W.  K.H. ,  and  Phillips,  M.  (1962)  Classical  Electricity  and  Mag¬ 

netism,  Addison-Wesley,  Reading,  2nd  Edition,  p.  407. 

12.  Mitchell,  A.  C.  G. ,  and  Zemansky,  M.W.  (1934)  Resonance  Radiation  and 

Excited  Atoms,  Cambridge  University  Press,  p.  116. 
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by  the  ground  atom  and  then  re-emitted.  This  process  is  quite  different  than  the 
scattering  process  assumed  in  the  present  paper.  Still  we  shall  compare  one  of 
the  consequences  of  these  assumptions  witii  that  of  the  present  paper  (scattering, 
not  absorption  and  re-emission)  as  given  by  Eq.  (8).  By  the  use  cf  detailed  bal¬ 
ance  and  relations  between  the  Einstein  A  and  B  coefficients  one  can  construct  a 
cross-section  E  considered  as  a  function  of  the  frequency  v.  1.  is  shown  that 


r  2 

g 


rdv  "  S  gTr  • 


(9) 


where  Xq  is  the  wavelength  of  the  resonant  line,  g,,  and  gj  are  the  degeneracies 
of  the  upper  and  lower  states,  and  T  is  the  lifetime  of  the  upper  stat"\  From 
Eq.  (8) 

Jud.  =  £  ^Edk  =  fg  \y\.  (10) 


But  K  -  (2ir/XQ).  From  the  Heisenberg  uncertainty  principle 

T  (AE)  —  "K  (n) 

where  AE  is  the  "width"  of  the  excited  state.  The  width  of  the  excited  state  is 
|  2y|  in  terms  of  wave  numbers.  But  AE  =-fic  |  2'y| .  Then  Eqs.  (10)  and  (11)  yield 
the  following: 

X  2 

f}T'du  ~~W  3  r'  ‘  <*2> 

Equation  (12)  agrees  with  the  phenomenological  result  Eq,  (9)1 

We  can  also  give  the  radiation  pattern  for  the  scattered  radiation.  We  shall 
give  the  case  where  the  incident  radiation  is  plane-polarized.  Let  i^,  i^  I3  be 
an  orthonormal  triad  of  unit  vectors.  We  take  the  incident  radiation  to  have  the 
form 

~in^~’  ^  =  Cii  sin  *  *“ct)  » 

H.n(x,t)  =  Ci_2  sin  k(i_3*x-ct)  .  (13) 


The  radiation  scattered  in  the  direction  described  by  the  unit  vector  %  is 
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by  the  ground  atom  and  then  re-emitted.  This  process  is  quite  different  than  the 
scattering  process  assumed  in  the  present  paper.  Still  we  shall  compare  one  of 
the  consequences  of  these  assumptions  with  that  of  the  present  paper  (scattering, 
not  absorption  and  re-emission)  as  giver  by  Eq.  (8).  By  the  use  cf  detailed  bal¬ 
ance  and  relations  between  the  Einstein  A  and  B  coefficients  one  can  construct  a 

f  i  o 

cross-section  T  considered  as  a  function  of  the  frequency  v.  is  shown  that 


hL 

*i T  ’ 


(9) 


where  \Q  is  the  wavelength  of  the  resonant  line,  gg  and  g^  are  the  degeneracies 

of  the  upper  and  lower  states,  and  r  is  the  lifetime  of  the  upper  stat"'.  From 

Eq.  (8) 

jldr  =  ±  jEdk  =  ff  \y  j  .  (10) 

But  K  =  (2 tt /XQ).  From  the  Heisenberg  uncertainty  principle 

r  (AE)  —  -h  (n) 

where  AE  is  the  "width"  of  the  excited  state.  The  width  of  the  excited  state  is 
|  2y|  in  terms  of  wave  numbers.  But  AE  =-fic  |  2'y| .  Then  Eqs.  (10)  and  (11)  yield 
the  following: 

A  2 

3  7  -  <121 

Equation  (12)  agrees  with  the  phenomenological  result  Eq.  (9)! 

We  can  also  give  the  radiation  pattern  for  the  scattered  radiation.  We  shall 
give  the  case  where  the  incident  radiation  is  plane-polarized.  Let  ^g  i.3  be 
an  orthonormal  triad  of  unit  vectors.  We  take  the  incident  radiation  to  have  the 
form 


Ein(x,  t)  =  Ci_j  sin  k(i_3  •  x-ct)  , 

Hin(x,t)  ■=  Ci,2  sin  k(i_3-x-ct)  .  (13) 


The  radiation  scattered  in  the  direction  described  by  the  unit  vector  T£  is 


10 


(14) 


Esc  (x.t)  =  [Uj  X  n)  X  nl  C(F(k)]1/2  cos[k(£.  x-  ct)  +  $(k)]  , 

Hgc  (x.t)  =  llfi-i  x  Dl  C(F(k))1/2  cos[k(2.-x  -  ct)+$(k>], 

where  4>(k)  is  defined  by 

sin*(k)  =-(F(k)]1/2  =  y  [(k-/£-6)2+ y2]"1/2  , 

cos  $(k)  =  (k-«-6)  [(k-/c-6)2  +  y2]'i/2.  (14a) 

At  resonance,  4>  =  -ff/2.  It  is  to  be  noted  that  |(i_^  X  q)  X  gj  =  li1  X  77 1  =  sin© 
where  0  is  the  angle  between  i^  and  £.  Hence,  the  maximum  radiation  is  in  the 
~2  *  -i-3  P^ane*  and  there  is  no  radiation  in  or  opposite  to  the  i^  direction. 

3.  EIGENFUNCTIONS  OF  THE  UNPERTURBED  HAMILTONIAN. 

EXACT  MATRIX  ELEMENTS 

For  the  sake  of  brevity  the  remainder  of  the  paper  will  be  written  as  a  direct 
5 

extension  of  earlier  work,  except  that  in  the  present  paper  the  helicity  will  be 

denoted  by  0  instead  of  A  as  in  5  to  prevent  confusion  with  the  wavelength  A  . 

5 

Equations  from  the  previous  paper  will  have  a  prime  attached  to  them. 

The  Hamiltonian  of  the  atom  interacting  with  the  radiation  field  is  the  usual 

or.e: 

H  =  HA  +  Hp  +  Hj  =  Hq  +  Hj,  (15) 

where  HA  is  the  usual  Hamiltonian  for  the  hydrogen  atom,  Hp  is  the  Hamiltonian 
for  the  photon  field  given  by  the  first  of  Eq.  (37'),  and  HT  is  the  interaction  of 
Eq.  (1).  As  is  .customary,  we  ignore  the  A  term.  Of  course,  Hq  =  HA  +  Hp  is 
the  unperturbed  Hamiltonian. 

The  space  of  wave  functions  is  sj.  med  by  the  eigenfunctions  of  Hq  .  Thes*- 
eigenfunctions  are  direct  products  of  atomic  eigenfunctions  and  either  vacuum 
states  of  the  field  or  of  n-photon  states  which  we  take  to  be  in  the  energy-  angular 
momentum  representation  as  discussed  in  Eqs.(4l')  -  (43').  Particular  eigen¬ 
functions  will  concern  us: 

(1)  The  state  in  which  the  atom  is  in  the  IS  state  and  the  photon  field  is  in  the 
vacuum  state  will  be  designated  by  1 1>.  The  state  in  which  the  atom  is  in  the  IS 
state  and  there  is  a  photon  in  the  state  whose  energy  is  Ep  =  fick,  whose  angular 
momentum  is  given  by  the  quantum  numbers  j,  m  and  whose  helicity  is  given  by 
/3  will  be  denoted  by  1 1,  E  ,  j ,  m ,  0> . 

r 
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(2)  The  analogous  kets  in  which  the  atom  is  in  the  2S  state  will  be  denoted 

by  |2>  and  |  2,  E  ,  j,m,/3>  respectively. 

P  3 

(3)  In  the  case  that  the  atom  is  in  a  2  P  state  with  the  magnetic  quantum 
number  M  (=  0 ,  ±  1)  and  the  photon  field  is  in  a  vacuum  state,  the  eigenket  will  be 
denoted  by  |  2 ,  M> .  When  a  photon  is  present  the  ket  will  be  written 
|2.M.Ep,j,m,0>. 

Thus  in  our  notation,  the  absence  or  presence  of  the  quantum  number  M  indi¬ 
cates  whether  the  atomic  state  is  an  S  state  or  P  state. 

These  three  sets  of  kets  are  orthogonal  to  all  other  eigenkets  of  Hq  .  They 
are  also  orthogonal  to  each  other.  Within  each  set  they  satisfy  the  orthonormality 
relations 


<11 1>  =  1 , 

<l,Ep.j,m^|l,Ep'.j'.m'.r>  =  Ep6(Ep-Ep')6..,  6^  6^,  , 

<l|  1  ,E  »j,m,0>  =  0, 

<2|2>  =  1, 

<2,Ep,j,m^|2,Ep',j',m',^>  =  Ep6  <Ep  -  Ep'>  6..,  6^,  6^,  , 

<  2 1  2  ,  E  ,  j ,  m ,  0  >  =  0, 

r 

<  2 ,  m|  2  ,  iVl'  >  =  6mm/, 

<2.M.Ep.j.m^|2,M'.Ep',j',m'.^>=6MM,  Ep 6  <Ep  -  Ep')  6.,  6^,  6^,  , 

<2,m|2,M',E  ,j,m,/3>  =  0.  (16) 

r 

We  shall  take  all  matrix  elements  of  Hj  to  be  zero  except  those  between  the 
eigenfunctions  introduced  above.  The  only  non-vanishing  matrix  elements  are 
easily  calculated  from  Eqs.  (2a')  and  (2b/)  when  one  uses  the  appropriate  radial 
wave  functions  for  the  hydrogen  atom.  The  non-vani3hing  matrix  elements  are  the 
following: 


:n,E  ,j,m,0  |Hj|2,M>  =  <2 ,  m|  Hjjn,  E  ,  j.,  m  ,/9>  * 


=  -i^(e2/a)(a/ir)1/9'6Mm6j>1Gn(ka). 


(17a) 


& 
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where  n  is  the  principal  quantum  number  1  or  2  and  the  functions  Gn(x)  are 
given  by 


G^x)  =  -  (2/3)1/2x[x2+ (3/2)2]'2  , 

G2(x)  =  -  (12f1/2x3[x2+  l]'3  .  (17b) 


Furthermore, 

<n|HI|2,M,Ep,j,m,jS>  =  <2 ,  M ,  Ep .  j ,  m.P  I  hJ  n>  * 

=  -(-l)mi|8(e2/a)(a/n)1/26M  _m6.ilGn(ka).  (17c) 


Finally, 

<2,M'|HI|2,M,Ep,j,m,)3>  =  <2 ,  M ,  Ep ,  j ,  m  ,  0 1  Hj|  2 ,  M'>  * 

■  .i-  i) 

X  6.  ,A(ka)  ,  (17d) 

j »  1 

where  the  function  A(x)  is  given  by 

A(x)  =  (2)"1/2x2[x2  +  l]"3  .  ( 17e) 

As  noted  previously  k  =  Ep/hc  . 


4.  EIGENFUNCTIONS  OF  THE  PERTURBED  HAMILTONIAN 
IN  THE  TWO  LEVEL  APPROXIMATION 


We  shall  denote  the  energy  of  the  ground  state  of  the  hydrogen  atom  by  Ej 
and  the  energy  of  the  first  excited  state  by  Eg.  Clearly  K  =  (Eg  -  Ej)/tic.  We 
shall  be  interested  in  particular  eigenfunctions  of  the  Hamiltonian  H  =  Hq  +  IT, . 
They  will  be  denot4d  by  |  1 ,  E  ,  j ,  m  ,0)  and  satisfy  the  set  of  integral  equations 
symbolically  given  by 


1  ,E  , j,m,0)  = 

r 


1  ,E  ,j,m,  j3>  +  y_  (Ej  +  E 


H0)  Hj|  1 


Ep,j,m,/3),  (18) 
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where  y_(x)  (which  is  not  to  be  confused  with  the  half-width  y  or  the  function  y(x) 
which  will  be  used  later  as  a  "half-wiath  function")]  ib  given  by 

y  (x)  =  ■-  -  iff6(x)  =  lim  —  ?? ■  ■  .  (19) 

x  t-*+0  e 

where  P  means  that  the  principal  value  of  the  integral  is  to  be  used  when  P/x  is 

used  as  a  factor  in  an  integrand.  The  eigenfunction  1 1 ,  E  ,  j ,  m ,  0)  is  an  outgoing 

13  P 

wave  function  in  the  sense  of  scattering  theory.  In  our  discussion  of  scattering, 

4  14 

we  shall,  however,  follow  more  closely  Friedrichs  and  Moses. 

It  is  readily  seen  that  j  1 ,  Ep ,  j ,  m ,  0)  satisfies  required  eigenvalue  equation 

(E-H0)  |l,Ep,j,m,0)  =  Hjll.E  ,j.m,0)  ,  (20) 

where 


E  =  E1  +  Ep. 


(20a) 


We  require  these  eigenfunctions  to  satisfy  Eq.  (18)  because  they  are  particu¬ 
larly  useful  in  solving  the  initial  value  problem  in  which  we  are  interested.  From 
Eq.  (18)  it  follows  that 


lim  exp  -jr  (Hn  -  EH 


|  1  ,Ep,  j,m,0) 


lim  exp  —  (Hq - E)t  |  1,E  ,j,m,0) 
t  ->  +0O  ** 


|  l,Ep,j,m,/3>  ,  (21) 

i  l,Ep,j,m,0> 

-2iri6(E-H0)Hj|l,Ep,j,m,|3).  (22) 


It  should  be  noted  that  the  right-hand  side  of  Eq.  (22)  is  an  eigenstate  of  Hq 
with  the  eigenvalue  E  .  It  can  be  shown  that  the  eigenfunctions  1 1 ,  Ep ,  j ,  m ,  0) 
satisfy  the  orthonormality  relations 


(l.E  ,j,m,p|l,E  j',  m',  0'}  = 


E  5  (E  -  E  ')  6..,  6  , 

P  P  P  31  mm 


(23) 


Let  us  require  that  the  solution  |  $  (t)  >  of  the  perturbed  time-dependent 
Schroedinger  equation  approach,  as  t-»-®°,  the  solution  of  the  unperturbed  time- 


13.  Lippmann,  B.A.,  and  Schwinger,  J.  (1950)  Variational  principles  for 

scattering,  I,  Phys.Rev.  79:469. 

- “ -  M\ 

14.  Moses,  H.  E.  (1955)  The  scattering  operator  and  the  adiabatic  theorem, 

Nuovo  Cimento  1:103. 

- m 
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dependent  equation  which  corresponds  to  the  state  in  which  the  atom  is  in  the 
graund  state  and  in  which  a  photon  is  in  an  energy,  angular  momentum  state.  The 
state  of  the  system  is  the  direct  product  of  the  ground  state  of  the  atom,  which  we 
take  to  be  normalized  to  unity  and  the  wave  function  of  the  photon  g(Ep  ,  j,  m  ,/3) 
whose  norm  is 


Z  f  lg<E  ,j,m,0)|2 
j ,  m ,  0  ' 'o  p 


Tr.en  for  t-*-°° 


|$>(t)>  =  exp  I  -  -jr  H0t  I  Z  \  |l  .E  ,j,m,0>g(E  .  j,m,/i>-g-2 
1  n  J  j.m.P  “o  v  p  p 


•  s  rii.Ep,j.m,«>0Xp[-i,E].Ep„] 


j.  m,0 

dE 

X  g(E  , j.m./J)  • 

*  n 


(24) 


For  t  finite  the  state  is  given  by 


l‘Ht>  >  =  exp  f-  xr  lit]  Z  {  |l,E  ,  j,m,j3)g(E  ,j,m,/S)-pr-2 
L  n  J  j.  m,  P  -!0  p  p  *> 

fcO  _  ,  _ 

ll.E  ,j.m,j8)  exp  j-  £  (E.+E  )t| 

j.m.p  -0  H  L  H  J 


dE 

X  g(E  .j.m.p) -g-E 

*  ■'n 


(25) 


For  t-*+« 


$(t)>  =  exp  k  I!ot]  !  f  I  l.E  ,j,m,P>g(E  ,j,m,P)-J^ 

1  J  (j.m.P  J0  P  P  LP 

i>00 

6(E-HJ  II.  |  l.E. j.m.p) 

0  0  I  P 

dFIpl 


X  g(Ep. j.m.p)  -g-Bj 


(26) 


)5 


where  E=E  +  .  as  before.  By  evaluating  <1 , E  ,  j , m , P  1 4»(t) >  using  Eq.(26), 

p  .  p 

it  is  seen  that  the  s:attei  ed  part  of  the  photon  wave  function  in  the  energy,  angular 
momentum  representation  is 


g..,  (E  =  -  2ffi  Z  <  1,  E,  j,  m,  ^  |  H.  1 1 ,  E,  j',  m'.  P') 

sc  p  j,  Ki.fr  H  H 

X  (Ep)~ 1  g  (Ep .  j',  m'.  P')  exp  [-  i  Ept]  . 


<  27) 


In  terms  of  the  linear  momentum  representation,  the  scattered  part  of  the  photon 
wave  function  is  given  by 


fsc(*-ftt) 


cl/2 

p 


ZZ  Yf‘^(0,<»g  (cp.j,  m,/3;t). 
j  m  J  SC 


(28) 


(See  Eq.  (31a7).]  The  angles  9  a.md  0  a  the  usual  angles  which  give  the  direc¬ 
tion  of  £  in  spherical  coordinates. 

At  finite  times  t  let  the  wave  function  h(£,/S;t)  be  defined  by 

1/2 

h(£,0;t)  -  f  ZZ  Y^(MXl,E.j,m,^k(t)>, 
p  j  m  J  p 

(Ep  =  cp).  (29) 

Let  £  =  (£/p)  and  define  v(£,£;t)  and  w(£,8;t)  by 

v(£,^t)  =  ^  rih(£,0;t)|2  p2 

w(£./3;t)  =  ^-  y  |h(£,^;t)|2p2  dp.  (30) 

Then  v(r},/3;t)dn  where  d.Q  =  sinO  d0d0 ,  represents  the  number  of  photons 
with  polarization  described  by  j3  passing,  per  unit  time,  through  the  solid  angle 
dfl.  Likewise  w(rj,^;t)  dft  gives  the  energy  of  radiation  per  unit  time  passing 
through  the  solid  angle  dfi  when  the  radiation  is  polarized  in  the  manner  given  by 
P .  The  atom  is  in  the  ground  state.  When  the  wave  function  of  the  incident  photon 
g(Ep,  j,m,/3)  is  chosen  to  correspond  in  the  limit  to  the  case  in  which  the  photon 
is  monochromatic  and  when  its  incident  path  is  a  prescribed  straight  line, 
v(j),)3;t)  and  w(£,/3;t)  are  independent  of  time.  One  can  then  use  these  quantities 
to  obtain  the  cross-sections.  Our  treatment  will  correspond  to  the  discussion  of 


16 


15  16 

time-pi  cportional  transition  probabilities,  *  which  in  turn  are  more  careful 

13 

treatments  of  the  theorv  discussed  by  Lippmann  and  Schwinger. 

li  is  seen  that  tc  solve  the  scatterinf  problem,  one  must  solve  Eq.  (18)  for 
the  eigenstates  1 1  ,  Ep,  j ,  m  ,/3).  Usually  these  eigenstates  are  obtained  by  pertur¬ 
bation  theory  o",  equivalently,  by  iteration.  On  the  right-hand  side  of  Eq.  (18), 
one  replace.*  1 1 ,  Ep ,  j ,  m ,  /J)  by  1 1 ,  Ep ,  j ,  m ,  /3>  and  continues  this  iterative  proc¬ 
ess  in  an  cbvious  way.  The  second  order  iteration  yields  results  which  are  of  the 
Kramers-Hoisenberg  dispersion  formula  type,  and  gives  infinite  results  when  the 

incoming  photon  has  its  energy  equal  to  the  difference  in  energies  of  two  atomic 

3 

states.  For  this  reason  Dirac  proposed  a  different  type  of  approximation  which 
is  more  accurate  near  resonance.  What  follows  is  our  version  of  the  Dirac  res¬ 
onance  theory.  Our  version  is  a  generalization  in  which  degeneracy  of  the  atomic 
state  is  taken  into  account.  Normally  one  would  need  to  use  analogues  of  "stabil¬ 
ized  eigenfunctions"  of  Hq.  But  because  we  are  working  in  an  angular  momentum 
representation  and  because  both  Hq  and  the  two  level  H  commute  with  the  angu¬ 
lar  momentum  (as  does  the  exact  H,  of  .ourse)  there  is  a  considerable  simplifica¬ 
tion.  Indeed,  the  use  of  the  correct  angular  momentum  representations  for 
photons  makes  this  comparatively  simple  calculation  possible. 

In  our  approximation  we  assume  that  the  only  non-zero  matrix  elements  of 
Rj  are  those  given  in  Eq.  (17).  These  are  all  the  one-photon  state  which  connect 
the  ground  state  and  the  first  excited  states  with  each  other.  The  effect  of  this 
assumption  is  that  we  have  replaced  the  interaction  Hj  by  another  interaction 
which  is  also  Hermitian.  This  interaction  enables  us  to  solve  Eq.  (18)  exactly. 
From  Eq.  (18) 


<l.Ep',j'.m',/J'|l.Ep.j.m./}>  =  V  <  V  Ep  >  fy  W  V 


+  y-  (Ed’Ed  1  5  <1'ED’’^'m'’^lHll2,M><2,M|l,Eri,  j,rn,j3).  (31) 


15.  Friedrichs,  K.O.  (1952)  Zur  asumptotischen  Beschreibun  von  Streuprozessen, 


/M 

16.  Moses,  H.  E.  (1953)  The  Scattering  Operator  in  Quantum  Mechanics,  Part  II: 
The  Scattering  Operator  Formalism  and  Other  Formalisms,  Research 
Report  No.  CX-13,  New  York  Univ. ,  Inst,  for  Mathematical  Sciences. 


Nachrichten  der  Akademie  der  Wissenshaften  in  Goettingen,  Mathematische 
Physikalische  Klasse,  Ila  Mathematish-Physikalish-Chemische  Abteilung  ” 
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On  using  Eq.  (17a) 

'  1 ,  Ep',  j',  m',  &'  1 1 ,  Ep ,  j ,  m ,  #  -  k6  (k-k')  6„/  6m_/6^/ 

-i(u3/J7)1/2a‘1i8,y  (k-k'jS./  ,0,  (k'a)<2 ,  m'l  1 , E^ ,  j ,  m ,0).  (32) 

In  Eq.  (32)  and  later  k=  Ep/?ic  and  k*  =  E^/fic. 

From  Eq.  (20)  we  have,  on  multiplying  through  by  <  2 ,  IVI  | 

p  00 

fe  -  (E,-E.)l  <2,m|i,E  .j,m,0)  =  ,  Z  \  <2,  M  |  IlJ  1 ,  E  '  j',  m',0'> 

I?  z  1  J  p  j’,  m',  0'  Jo  P 

dE  ' 

X  <1,E  '  j'.m'.fl'Il.E,  ,j.m,0>  — f 

P  P  E 

P 

ico 

<  2 .  M I  H.  I  2 ,  E  7,  j',  m'(  0'> 

1  D  J 

j»*«  0 

dE  ' 

X  <2,E  ',j',m'.0'|l,En,j.m,0)-_E. 

P  P  Ep 

,■*00 

+  Z  \  < 2 ,  M  |  H.  1 2 ,  iVl',  E  /,  m',  0'  > 

j',m',0'  *>0  1  p 

M' 

.  dE' 

X  <2,M',E  '  j',m',0'|  l,E.j,m,0) — p  .  (33) 

p  p  V 

On  using  the  expressions  Eq.  (17)  for  the  matrix  elements 

(k-K)  <2,m|  1  ,E.j,m,0)  =  i(a3/7/)1/2a-1  Z  0' f  G.(k'a) 

p  0'  ''0  1 

X  <  1 ,  Ep' ,  1  ,  M  ,  0'  |  1  »  Ep ,  j ,  m ,  0)  j~r 

+  i(o3/ir)1^2a'1Z/3,J  G2(k,a)<2,Ep/.l,iVI,/S'|l,Ep,j.m,^)^- 

+  (-t)M^1(«3M>1/2a-1  I  f1  1  1  \ 

m’,  fi1  \m'  -M  M-m'/ 

X  f  A(k/a)<2,M-m'.E  l.m'.^'ll  ,E  ,j,m,(3)  ^  ,  (34) 

PH 


where  we  have  used  k 


e2*ei 


18 


From  Eq.  (18)  one  obtains  in  a  similar  fashion  the  following: 


<2,Ep\j',m',/J'|l  E  .j.m.JJ)  = 

-  i6.,  ,/3'(a3/ir)1/2a'1-y  (k-k'-K)G9  (k'a)<2 ,  m' | 1  ,E  .  j,m,/3),  (35) 

<2 ,  M', E  ',  j',  m'./3'l  1 ,  E  ,  j ,  m  ,  fi)  = 

P  P  * 

6,(1ta3W1/2a‘1y.  (k-k'-/c)  Ji^(-l)m/6m,  _M,G1(k,a) 

X  <l|  1  ,Ep, j,m,0) 

+  i0'(- l)m'  6m,  _M,  G2 (k'a)< 2|  1 ,  Ep ,  j ,  m ,  fi) 

-(  /  ,  1  j  (-l)M,+  m/ A(k/a)<2,m,4  M'|  l,En,j,m,p)j,  (36) 

\m  -(m'+M')  MV  P  ) 

<l|l,E  ,Jm,0)  =  -i(a3/ir)1/2a'1y.(k) 

X  Z  fl'(-l)mT  G1(k'a)<2,-m'E  '  1  .m'.P'l  1  ,E.  j ,  m,/H$r,  (37) 

m  ,0'  J0  P  p  k 

<2|l.E  .j.m./J)  -i(a3/ir)  i/2a"1y_  (k-K) 

/  pOO 

X  L  \  G0  (k,a)<2 ,  -m' ,  )3'  1 1  ,  E  ,  j  ,  m  ,  /3) 

m ,  j3'  2  P 

X  •  (38) 

Equations  (32),  (34),  (35),  (36),  (37),  and  (38)  are  a  set  of  equations  for  the 

quantities  <  1 ,  Ep',  j',  m',0'1 1  ,Ep,  j ,  m , 0);  <2 ,  m|  1  ,Ep  ,j ,  m , jj); 

<2  ,  Ep',  j',  m',/}'  1 1 ,  Ep,  j  ,  m  ,0);  <2,M',  j',  m',/3'|  1 ,  Ep ,  j  ,  m  ,  0); 

<  l|  1 ,  E  ,  j ,  m ,  /3);  and  <2|  1 ,  E  ,  j ,  m ,  0).  It  is  to  be  noted  that  these  are  the 
*  I 

only  components  of  |  1 ,  E  ,  j ,  m ,  0)  in  the  Hq- representation  which  have  the  pos¬ 
sibility  of  being  non-zero. 

We  shall  now  solve  the  set  of  equations.  On  substituting  Eq.  (36)  into  Eqs. 

(37)  and  (38),  one  obtains  a  pair  of  homogeneous  linear  equations  for 
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<l|  l,Ep,  j,m,0)  and  <2|  1 ,  Ep,  j ,  m ,  fi).  Since  the  determinant  of  the  coeffi¬ 
cients  is  not  zero  it  follows  that 

<l!  l,Ep,  j,m,/3)  =  <2|l,Ep>j.m,/3)  =  0.  (39) 

Equation  (36)  now  simplifies  considerably.  On  substituting  Eqs.  (36),  (32), 

and  (35)  into  Eq.  (34),  one  obtains  an  equation  for  <2 ,  m|  1 ,  Ep,  j,  m,/?)  which 

one  can  solve  easily.  On  substituting  <  2 ,  M 1 1 ,  Ep ,  j ,  m ,  /3)  so  obtained  into 

Eqs.  (32),  (35),  and  (36)  one  finds  <  1 ,  Ep',  j',  m' ;fi'  1 1 ,  Ep,  jm ,/?), 

<2,E  ',j',  m'.jS'l  1 ,  E  .  j ,  m ,  j3),  and  <2 ,  M',  E  j',  m'/S'  1 1  ,E  ,j,m,/3)  respec- 
P  P  P  P  j 

tively.  In  carrying  out  the  substitutions,  one  should  note  that  summations  over  p 

yield  factors  of  2.  Furthermore,  one  uses 


if1  1  1  ' 

M'  \m'  -M  M-m' 


for  M  =  0 ,±  1 . 

In  order  to  give  the  solution,  it  is  convenient  to  make  some  definitions.  Let 
the  functions  L(x)  be  defined  for  i  =  2,  3  by 


_  oo  w 

L(x)  =  J  jG.(|)j2y_  (x-£)  for  i=  1,2. 
I3(x)  =  J  [A(?)]2y_(x-?)  f-. 


The  functions  Ij(x)  are  given  in  Appendix  I.  Indeed  the  imaginary  pares  can 
be  obtained  immediately  using  the  5-function  in  Eq.(19). 

Let  us  fuither  define 

2a2  2ot2 

«j(k'  =  —  Re  Ij(ka),  $2(k)  =  rra  Re  Igfka-Ka), 


63(k)  =  3ifa  Rel3(ka-Ka). 


Let  H(k)  =  1  for  kSO  and  H(k)  =  0  for  k<0.  Then  also 


2a3  -  .  „  ,  2a3 

-z—  Im  I.(ka)  = - 5 

ffa  1  ka2 


ygfk)  =  Im  Igfka-Ka) 


(k-K^a 


2  ^Gg(ka-Ka)J2  H(k-K), 
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Ort3  «-y3 

y..(k)  =  ~rr  Im  IJka-fCa)  =  - ■= 

J  *  Stk-'Oa'3 


j^A(ka-Ka)j2  H(k-/c) 


3  3 

6<k)  =  Z  6.(k)  .  y(k)  =  Z  y.(k) . 

i=  1  1  i=l  1 

[6(k)  of  Eq.  (45)  should  not  be  confused  with  usual  6-function.] 


<l,Ep,,j/,m,./S'|l.Ep.j,m.fl)  =  kfi(k-k')  6.  ,  6mm,  6^, 

+  2?  6j,l  6;,j'  6m,m/  y-(k-k,)  k-  k-600  -  iy(k) 

X  [kk'^My^k')]1'2  . 

.  ^[-(k/2tr)y.(k>]1/2 

<2,M|l.Ep.i.m.»  .  «M>n,  ■ 

<2 '  Ep  '  ^  1 1 '  Ep 1  i ' m  =  2i7  "  j ,  1  ^j  j'  PP>7_  (k-k'- 


[kkVj(k)  y2(k'+>c)]1/2 

k-K-6(k)  -  iy(k)  ’ 


<2,M,.Ep,.j,.m,.^|l.Ep,j.m,/3)  = 


'  /  I  11 


J*x  -m  mV 

[skk'y^k)  y3(k'+K)]1/2 

X  k-K-6(k)  -  iy(k)  '  1 

A  quantity  which  is  important  for  finding  the  scattered  portion  of  the  wave 
function  [see  Eq.  (27)]  is  given  by 

-2iri(Ep)'1  <l.Ep.j.m,j3|H1|lfEp.j,.m,(/3')  = 

-»i(ttck)_1  Z  <  1 ,  E  ,  j  ,  m  ,  j3 1  Hj I  2  ,  M ><  2 ,  M 1 1 ,  E  j' ,  m',  0')  = 
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(50) 


6j,l  6j,j'  6m,m'  ^’rl(k)  k- /e-6(k)  -  iy(k)  * 


It  is  to  be  noted  that  angular  momentum  is  conserved  as  indicated  by  the  pres¬ 
ence  of  factors  as  6.  ./  6 _ / .  Eq.  (39)  is  another  consequence  of  the  conserva- 

J.J  m,m 

tion  of  angular  momentum. 


5.  THE  LINE  BREADTH  AND  THE  SHIFT  IN  POSITION  OF  RESONANCE 

It  is  obvious  that  the  square  of  the  absolute  value  of  k-*c-6(k)  -  iy(k)  -1,  that 
l  2  2|-1 

is  j  [k-/c-6(kl]  +  [y(kl]  j  ,  is  tl  •>  resonance  denominator  in  the  scattering  cal¬ 

culations.  Since  we  are  considering  values  of  k  near  k  ,  and  since  6(k)  and  y(k) 
are  slowly  varying  functions  of  k.  we  may  replace  k  by  K  in  these  functions  as 
is  customary.  We  define  6  and  y  by 

6  =  6 (K),  y=y(K)  s  y^K).  (51) 

The  res.onance  denominator  then  becomes  [(k-K-6)2  +  y2]"1.  The  values  6 
and  y  given  by  Eq.  (5a)  are  those  obtained  from  Eq.  (51).  It  is  seen  that  while  the 
matrix  elements  corresponding  to  the  "forbidden"  transitions  contribute  nothing 
to  y,  they  give  a  substantial  contribution  to  6.  The  value  of  6  in  6,7,  which  ig¬ 
nores  the  effect  of  the  forbidden  transitions  is  the  principal  contribution  of  6j 
which  is  seen  to  be  6^0).  The  value  of  y  of  6  is  in  error  and  should  have  the 
value  given  in  Eq.  (5a). 


6.  CALCULATION  OF  THE  CROSS-SECTIONS 


We  shf!1  now  derive  Eqs.  (4)  and  (6)  for  <7j(0)  and  crg^  respectively.  Our 
procedure  will  be  to  obtain  the  transitions  per  unit  time  w(r?  ,(3;t)  of  Eq.  (30), 
taking  as  the  initial  photon  state  g(E  ,  j,m,/S)  that  one  which  gives  rise  to  a  cir- 

r 

cularly  polarized  electromagnetic  wave  of  a  given  wavelength,  travelling  along  the 
positive  z-axis.  One  then  divides  w(?7 ,  /3;  t'  by  the  flux  of  energy  of  the  incoming 
radiation,  that  is,  the  Poynting  vector. 

We  shall  now  proceed.  It  is  convenient  to  specify  the  initial  photon  state  in 
the  linear  momentum  representation  instead  of  the  angular  momentum  represen¬ 
tation.  Let  the  initial  state  in  the  linear  momentum  representation  be  denoted  by 
f(g,)3).  Then  from  Eq.  (3 lb') 


g(Ep,j,m,/3)  = 


JT2 


nir 

1  dM 


tie  sin  6 
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p  =  (Ep/c)  (sin  8  cos  0,  sin  6  sin  0,  cos  0) . 

On  using  Eqs.  (29),  (46),  and  (52)  we  obtain  the  wave  function  of  the  photon  at 
finite  times  in  the  linear  momentum  representation  as 

Wfi.fct)  =  e'icktf(E,/3>  +  I  £  C  %  Y™''0(e,0> 

*  m,p  ^  p 

X  (0',  0')/3/3,G(k',k)  e”ick/^y_(k,,-k)  ffp'.jS') , 

where  0/,0/  are  the  polar  angles  of  £',G(k',k)  is  given  by 
[kk'y(k)y  1(k,)ll/2 

G(k'k)  =  2 it  k'-K-6(k')-iy(k;)  '  ( 


and  where  here  and  later 


E  =  cp  =  cfik,  E  '  =  cp'  =  cfik',  E  "  =  cfik"  .  (53 

P  P  P 

In  deriving  Eq,  (53),  we  have  used  the  completeness  relation  for  Y”1,m  (0,0) 
Eq.(A12).17 

For  the  wave  function  f  (£,/))  let  us  take 

f(E,/3)  =  6^  gf(p)  (5 

where  f(p)  is  a  real  function  of  £  which  has  a  very  sharp  peak  at  p  =  p^,  where 
£o=  fiK(0,0, 1),  (54 


The  incoming  electromagnetic  field  corresponding  to  this  wave  function  can  be 
obtained  from  Eq.  (34').  Because  of  the  sharpness  of  the  peak  in  the  wave  function 


E(x,  t)  =  - 


[&5<£0>  ciK,z-“’  -tefaj 


17.  Moses,  H.  E. ,  and  Quesada,  A.  F.  (1973)  The  expansion  of  physical  quan¬ 
tities  in  terms  of  the  irreducible  representations  of  the  sr-ale-Euclidean 
group  and  application  to  the  construction  of  scale-invariant  correlation 
functions.  II.  Three-dimensional  problems.  Generalization  of  the 
Helmholtz  vector  decomposition  theorem.  Arch.  Rati.  Mech.  50:194. 
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rrc  n  A  ,  iK(z-ct)  n*,  ,  .-iK(z-ct) 

S~  '  W72  %(£0  e  +  %(£o  e 


=  J  fte)d£- 


In  terms  of  the  Cartesian  components 


E  (x ,  t)  =  H  (x ,  t)  =  —rp) - 37o  sinfK(z-et)] 

AB 

E  (x , t)  =  H  (x,t)  =  — yj7) - -To  cos  [K(z-ct)] 

>  ~  x~  2r  ^  vYi  * 


E  (x.t)  =  H  (x.t)  =  0. 


(56) 


Only  the  z-component  of  the  Poynting  vector  is  not  zero.  The  z-component, 
denoted  by  S,  is  given  by 


S  =  E.  H  -  E  H 


cA 


x  y  y  x  4n 


(2ff«) 


3  * 


(57) 


We  shall  calculate  w(rj,/3;t)  for  T)  which  does  not  coincide  with  the  z-axis. 

For  this  reason,  and  because  f(p)  is  so  sharply  peaked  at  p  =  prt,  the  first  term 

—  ~  •‘-0 
on  the  right  of  Eq.  (53)  is  zero.  Thus  one  has 

p  m  m  "  p  "  p 

X  Y ™,,^*(0,<J)  Y G(k',k)  Y™  ’ho",#)  G*(k\k) 

x  (k'-W,)e-ic(k,“1^)tx.(k/-k)y*(k'-k)f(p/)f(£)  , 


(58) 


where  p7,  O'.ih',  and  p”,  0",  <j>"  are  the  spherical  coordinates  of  p  and  £  re¬ 
spectively.  In  Eq.  (5  we  make  the  substitution 

(k'-k')y.(k'-k)yV-k)  =  [(k'-k)-(k"-k)]  yjk'-k)  y_V-k) 

=  y*  (k"-k)  -  y_(k'-k) 


(59) 


and  obtain  two  terms  on  the  right  of  Eq.  (58)  On  interchanging  the  primed  and 
double-primed  dummy  variables  in  one  of  the  terms 
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&  lh<p.^;  t)i7'  -  -  T  to  Z  2  J  %  ^  »“'•*(«,<» 

p  m  m  “  p  “  p 

X  Ym'‘^*(9,  <j>)  Ym,’P*(e',6')G(k'.k)  Ym’^(0',0')G*(k',k) 


X  e'ic(k,*k*)ty_(k,-k)f(£  )  f(£) 


(60) 


We  now  use  the  fact  that  f(p)  is  sharply  peaked  at  g=  g^  to  make  the  following 
approximations  under  the  integrand  which  become  increasingly  good  with  increas¬ 
ing  sharpness  of  the  peak. 

-ic(k'-k')t  -ic(K-K)t  , 
e  ~e  =  l  , 

Ym^AV1<)')~Y1m,’^(0^')  =  (3/4tr)1/26m, 

Y™’P  (6",6’)~  Y™’$  (0,0")  =  (3/4fr)1/26m.  £  , 

p'~  fiK  ,  p'-hK  , 

G(k'.k) -G(K.k),  G(k',k)  ~G(K,k),  y.(k'-k)  ~  y_(K-k) .  (61> 

Then,  on  using  Im  y_(x)  =  -ff6(x) 

2  ^ 

|h(p.0;t|2  =  |g(K,  k)  |2  I ’  ^(9,  6)  I  2  6(K-k) .  (62) 

~  2p"(hK)4  1 

Finally,  from  Eq.  (30)  and  the  fact  that  6(K-k)  = 

w(g,/3;t)  =  j*  |h(g,)3;t)|2  p2dp 

3cA2  (y1(K)]2 

K28jt2«3  (K-tc-6(K)]2  +  |y(K)]2 

where  9. 6  are,  of  course,  the  polar  angles  of  rj . 

Now  from  Appendix  D  of  Ref.  17 

Y^(9,0)  =  (3/  16ir)1/2  (1  +  cos  9)  , 


6(p-  K) 


A  . 

I  Y^'*  (0.5)1.  (63) 


fa 
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(64) 


Y =  <3/16ir)1/2e2i^(l-cos  0)  . 

/*■*!- 

Now  the  cross  section  for  the  case  that  the  circular  polarization  of  the  scat¬ 
tered  photon  is  the  same  as  that  of  the  incident  photon  is  given  by 

9  Iri<K>]2  4 

oA6)  =  w(t) ,0;t)/S  =  -2-  - L— - n-  cos  10 /2)  .  (65) 

1  ~  4K^  [K-<C-fi(K)r+ [y(K)r 


Similarly,  the  cross  section  for  the  case  that  the  circular  polarization  of  the 
scattered  photon  is  opposite  to  that  of  the  incident  photon  is 


a =  w(tj  ,  -/5;  t)/S  = 


(Y^K)]2 


4K 


{K-k-600]2  +  [y(K)]2 


sin’(0/2)  . 


(66) 


The  expressions  Eqs.  (4)  and  (6)  are  obtained  from  Eqs.(65)  and  (66)  re- 

2 

spectively  byreplacing  K  in  the  functions  K  ,  6(K),y(K),  and  y^K)  by  >c ,  and  by 
using  Eq.  (51), 


7.  CALCULATION  OF  THE  SCATTERED  ELECTROMAGNETIC  FIELD 

We  shall  now  derive  Eq.  (14)  for  the  scattered  electromagnetic  field. 

As  before  the  wave  function  of  the  incident  photon  in  the  linear  momentum, 
helicity  representation  will  be  denoted  by  the  complex  function  f<£,0).  On  using 
Eqs.  (27),  (28),  (50),  and  (52),  the  wave  function  of  the  scattered  photon  in  the 
linear  momentum  representation  is 

.  .  -  .  *y><k)  -ickt 

fsc(E'*:t)  =  1  k-x-soo- iyfi<7  e 

xIZj1  I*2*  d*'  C*  sin  9'  yT 

m  ff  J0  J0 

X  f<£',0')d (67) 

where  9,  $  and  9',$'  are  the  polar  angles  of  £  and  £;  respectively  and  |£!  * 

I  £*  I  =  P  =  «k. 
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The  incident  wave  function  will  be  taken  as 


f(p,j3)  =  f(p) 


and  is  independent  of  /3 .  The  function  f(p)  is  real  and  has  a  sharp  peak  at 
P  ~Pq»  where  p^  =  hK(0, 0, 1),  as  before.  It  is  readily  shown  that  the  incident 
electromagnetic  field  is  given  by  Eq.  (13)  where 


J  f(p)  djs . 


The  proof  is  close  to  that  for  Eqs.(55)  and  (56). 

The  electromagnetic  field  of  the  scattered  radiation  is  obtained  by  substituting 
f  (£,/?; t)  for  g(p,X)  er  *(i/fi)cpt]  in  Eq.(24').  The  scattered  electromagnetic 
field  can  be  written 


E(x,t)  =  ^  E(x,t;£)dn, 
H(x.t)  =  ^  H(x.t;2>cK2, 


o 

where  in  Eq.  (24')  we  have  written  dp  =  p  dpdft,  with  dfl=  sin  9  dO  d0  being  the 
element  of  solid  angle  and  £=■  £/ p  as  before.  The  vectors  Elx.tjj.)  and  H(x,t;£) 
are  the  components  of  the  electrofield  propagating  in  the  direction  £.  The  quanti¬ 
ties  E(x,  t;Tj)dQ  and  H(x,  t;£)dn  are  portions  of  the  electromagnetic  field  which 
would  be  intercepted  by  an  antenna  of  aperture  dfl.  In  fact,  E(x,  t ;£>  and  H(x,  t;Tj) 
are  Esc(x,t)  and  Hsc<x,  t)  of  Eq.  (14).  The  evaluation  of  E(x,t ;£)  and  H(x,  t;£)  in¬ 
volves  an  integration  over  the  variable  p.  Because  of  the  sharpness  of  f(jg)  at  p 
near  p^.  the  following  approximations  are  made: 


aik(£-x-ct)  iK(£«  x- ct) 


y\ 

k-K-C(k)  -  iy(k) 


n(K) 

~  K-/C-fi(K)-iy(K) 

=  -(F(K)l1/2ei${K) 


7 

K-fC-6-iy 


In  Eq.  (71)  we  have  used  the  fact  that  K  is  near  K  .  Eq.  (14a)  and  (51)  have 
also  been  used. 

The  sharpness  of  f(£)  also  permits  the  approximation  given  by  the  second  of 
Eq.  (61)  where  0  is  replaced  by  /3',  Finally,  one  has 
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^ - 


E(x , t; tj)  =  -(3/8ff)1/2C[F(K)]1/2  Re  cl[K(a*  ct  +  $(K>1 


X 


Z  Z 

fi=±i  p=±i  1 


(72) 


But  by  explicit  calculation 

Z  Z  ^'QJa.)  Yf,j3(0,0)  =  -  (3/2tr)1/2  [(i.Xn)  X  n].  (73) 

0=±1  /3'=±1  1  1 

The  first  of  Eq.  (14)  then  follows.  The  second  of  Eq.  (14)  is  obtained  in  a 
similar  fashion. 

8.  COMPLETENESS  OF  THE  EIGENSTATES. 

HIGHER  ORDER  PERTURBATION  THEORY 

By  eliminating  all  but  the  two-level  matrix  elements  of  the  interaction,  we  have 
constructed  an  approximate  Hamiltonian  H*.  One  can  ask  the  question  whether  the 
eigenfunctions  of  Section  4  are  a  complete  set  of  eigenstates.  We  believe  that  this 
question  can  be  answered,  and  missing  eigenstates  can  be  found  using  techniques  of 
Reference  4.  When  we  have  such  a  complete  set,  we  may  get  a  better  approximation 
by  considering  the  exact  Hamiltonian  H  to  be  a  sum  of  the  approximate  Hamiltonian  H', 
whose  eigenstates  are  known,  and  a  perturbation  consisting  of  the  interaction  which 
has  been  ignored  so  far.  One  can  set  up  integral  equations  for  the  "outgoing"  eigen¬ 
functions  of  H  in  terms  of  the  "outgoing"  eigenfunctions  of  H/.*8  This  integral 
equation  can  be  solved,  formally  at  least,  in  terms  of  a  Born  expansion  about  the 
"outgoing"  eigenfunctions  of  H/.  We  believe  that  the  expressions  for  these  Born 
approximation  terms  will  converge  if  the  exact  matrix  elements  of  the  electro¬ 
magnetic  interaction  are  used  to  prevent  ultraviolet  catastrophes.  One  of  the 
principal  causes  of  divergence  in  the  usual  Born  expansion  of  the  problem  is  that 
it  is  not  recognized  that  Hq  has  point  eignevalues  embedded  in  the  continuum 
which  the  interaction  may  cause  to  disappear.  Even  when  one  uses  H',  the  Born 
expansion  about  Hn  would  yield  divergences.  This  matter  is  discussed  in  great 

u  4 

detail  for  a  simpler  model.  Thus  a  more  careful  study  of  the  nature  of  the 
changes  of  the  spectrum  is  needed  to  prevent  some  of  the  divergences. 


18.  Gell-man,  M. ,  and  Goldberger,  M.  L.  (1953)  The  formal  theory  of  scattering, 
Phys.  Rev.  91:398. 
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Appendix  A 

Evaluation  of  the  Integrals  l;(x) 


We  shall  first  give  the  integrals  I.(x)  and  then  sketch  the  derivation  of  Ig(x) 
which  is  the  most  complicated  of  them,  the  derivation  of  I j  (x)  and  I^(x)  being  a 
simpler  version  of  that  for  Igix). 

We  need  give  only  the  real  parts  of  L(x),  since  the  imaginary  parts  are  given 
in  Eq.  (44). 
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It  should  be  noted  that  only  the  first  terms  on  the  right  of  Eqs.  (A-2)  and  (A-3) 
contribute  to  6,  since  x  =  0  in  the  computation. 

To  derive  the  expressions  for  I.(x)  the  following  integrals  are  used: 
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where  z  is  a  complex  number. 
Now,  from  Eq.  (41) 
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where  z  =  x+  i£.  After  evaluating  the  integral  we  take  6  -*0 , 
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